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Abstract 


This  paper  describes  a  technique  to  determine  interlaminar  tension  strength,  gqc,  of  a 
fiber  reinforced  composite  material  using  a  curved  beam.  The  specimen  was  a  unidirectional 
cunred  beam,  bent  90°,  with  straight  arms.  Attached  to  each  arm  was  a  hinged  loading 
mechanism  which  was  held  by  the  grips  of  a  tensile  testing  machine.  Geometry  effects  of  the 
specimen,  including  the  effects  of  loading  arm  length,  inner  radius,  thickness,  and  width,  were 
studied.  The  data  sets  fell  into  two  categories:  low  strength  corresponding  to  a  macroscopic 
flaw  related  failure  and  high  strength  corresponding  to  a  microscopic  flaw  related  failure.  From 
the  data  available,  the  loading  arm  length  had  no  effect  on  g3c.  The  inner  radius  was  not 
expected  to  have  a  significant  effect  on  03c,  but  this  conclusion  could  not  be  confirmed 
because  of  differences  in  laminate  quality  for  each  curve  geometry.  The  thicker  specimens  had 
the  lowest  value  of  030  because  of  poor  laminate  quality.  Width  was  found  to  affect  the  value 
of  030  only  slightly.  The  wider  specimens  generally  had  a  slightly  lower  strength  since  more 
material  was  under  high  stress  and,  hence,  had  a  larger  probability  of  containing  a  significant 
flaw. 


Introduction 

Because  of  low  interlaminar  strengths  of  laminated  composites,  interlaminar  failures  can 
be  a  predominant  failure  mode.  Interlaminar  failure  may  be  predicted  using  strain  energy 
release  rate  and  interlaminar  fracture  toughness  data  if  stress  singularities  are  present,  or,  if  the 
stresses  are  finite,  stress  and  interlaminar  strength  data  may  be  used.  In  the  latter  case, 
transverse  width  strength,  02c,  determined  from  fiat  90°  specimens  is  often  used  to  represent 
interlaminar  tension  strength,  However,  if  the  actual  interlaminar  tension  strength  is 
significantly  cfifferent  than  the  transverse  width  strength,  predictions  will  be  incorrect. 
Consequently,  a  method  is  needed  to  measure  the  interlaminar  tension  strength. 


Several  attempts  have  been  made  to  design  an  interlaminar  tension  strength 
specimen.  In  Ref.  1 ,  sixteen  unidirectional  specimens  of  24  plies  each  were  adhesively 
bonded  together  to  create  a  384-ply  bonded  laminate.  A  radius  was  machined  into  the 
specimen  to  give  a  minimum  area  in  the  center.  Aluminum  shanks  were  also  bonded  to  the 
ends  so  that  a  tension  load  couid  be  applied.  For  an  XAS/914C  carbon/epoxy  specimen,  an 
average  030  of  75.0  MPa  was  measured.  This  compared  with  an  average  020  of  83.0  MPa 
determined  from  flat  90°  specimens.  However,  failure  in  the  interlaminar  tension  specimens 
was  often  dose  to  a  bond  line.  In  Ref.  2,  a  similar  method  was  attempted  using  thick  (50  and 
100  plies)  unidirectional  and  cross-ply  laminates.  A  radius  was  machined  into  the  spedmens 
to  give  a  minimum  area  at  the  center,  and  aluminum  shanks  were  bonded  to  the  ends  for  load 
application.  For  an  AS4/3501-6  carbon/epoxy  specimen,  an  average  03c  of  43.0  MPa  was 
measured.  This  compared  to  an  average  02c  value  of  57.6  MPa  from  flat  [90]ioo  spedmens. 
However,  the  manufacture  of  thick  laminates  and  the  machining  were  considered  by  the 
authors  to  be  disadvantages  to  this  type  of  specimen. 

in  Ref.  3,  unidirectional  curved  laminates  were  used  to  determine  030  data  for 
AS4/3501  -6.  The  curved  laminates  were  essentially  "L"  shaped  with  the  interlaminar  tension 
failure  occurring  around  the  angle.  A  spedal  loading  fixture  had  to  be  manufactured  to  apply  a 
load  normal  to  the  loading  arm  to  open  the  specimen.  An  aveiage  03c  of  45.0  MPa  was 
measured  with  this  type  of  spedmen.  This  compared  to  an  average  020  value  of  65.4  MPa 
obtained  from  [90]24  flat  laminates.  In  Ref.  4,  a  semi-drcuiar  curved  beam  spedmen  was  used 
to  determine  03c  of  G40-600/5245C  carbon/epoxy.  These  laminates  were  of  poor  quality 
and  had  a  high  void  content  and  variations  in  thickness.  The  data  from  the  semidrcular  beam 
were  grouped  into  two  sets:  one  with  an  average  03c  of  58.0  MPa  and  the  other  with  an 
average  strength  of  32.5  MPa.  The  weaker  data  set  was  considered  to  be  a  flaw  related 
failure  whereas  the  stronger  data  set  was  proposed  to  be  a  true  material  property.  Also  in 
Ref.  4,  an  elliptical  curved  beam  spedmen  was  used  to  determine  030  of  T300/934 
carbon/epoxy.  The  manufacturing  method  was  improved,  and  a  very  high  quality  part  was 
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produced  with  no  discernable  flaws.  For  these  specimens,  a  mean  strength  of  107.1  MPa 
was  measured  which  was  almost  twice  the  in-piane  transverse  strength.  A  review  of  curved 
beam  testing  and  interlaminar  tensile  failure  was  given  in  Ref.  5. 

This  paper  describes  an  aitemate  configuration  of  the  cun/ed  beam  specimen  which 
offers  many  advantages  over  the  other  methods  for  determining  the  interlaminar  tension 
strength.  The  specimen  geometry  is  much  easier  to  manufacture  than  the  other  specimens. 
The  geometry  resembies  many  parts  currently  being  used  on  actual  aircraft  such  as  webs  and 
spars.  Also,  the  test  can  be  performed  in  a  simple  tensile  testing  machine  without  a  spedai 
loading  fixture.  The  test  section  of  the  curved  beam  is  more  representative  of  a  composite 
structure  than  a  flat  tensile  coupon  and,  hence,  may  give  more  representative  strengths  to  use 
in  faiiure  predictions. 

Specimen  Configuration  and  Test  Procedure 

The  curved  beam  specimen  is  shown  schematically  in  Fig.  1 .  Various  widths  w, 
thicknesses  t,  inner  radii  ri,  and  ioading  arm  iengths  L  were  used  as  shown  in  Table  1 .  The 
material  used  was  AS4/3501-6  carbon/epoxy.  The  material  properties  were  taken  from  Ref. 
3  and  were: 


Ee  =140GPa.  Er=11.0GPa,  Gre  =5.84GPa,  VtQ  =0.0237 


where  the  subscript  0  coincides  with  the  fiber  direction.  The  laminates  were  laid  up  in  300  mm 
wide  strips  over  the  comer  of  a  solid  aluminum  block.  The  comer  of  the  block  had  the 
appropriate  radius  to  form  the  inner  radius  indicated  in  Table  1 .  Each  of  the  straight  arms  was 
80  mm  long.  All  the  laminates  were  unidirectional  with  the  fibers  running  in  the  direction  shown 
in  Fig.  1 .  The  laminates  were  unidirectional  to  prevent  matrix  cracks  and  edge  stresses  from 
initiating  delaminations  [3].  The  laminates  were  ojred  in  an  autoclave  according  to  the  material 
manufacturer's  instructions.  A  thermal  blanket  was  also  used  on  the  48-ply  specimens  to 
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ensure  uniform  heating.  The  specimens  were  then  machined  to  the  specified  width  from  the 
300  mm  panel.  The  exact  width  and  thickness  of  each  specimen  were  then  measured  with  a 
vernier  caliper.  All  spedmens  were  dried  prior  to  testing  using  the  following  cycle:  1  hour  at 
95°C,  1  hour  at  1 10°C,  16  hours  at  125‘’C,  and  1  hour  at  150°C. 

initially,  loads  were  applied  via  aluminum  hinges  damped  to  the  loading  arms. 
Adhesively  bonding  the  hinges  to  the  surface  was  attempted,  but  the  adhesive  bond 
invariably  failed  prior  to  interlaminar  tension  failure.  In  some  configurations  of  the  beam,  the 
failure  loads  were  suffidently  high  to  bend  or  even  break  the  aluminum  hinges.  Hence,  a  steel 
loading  fixture  was  manufactured  to  replace  the  hinges  and  is  shown  in  Figs.  1  and  2.  The 
method  of  load  application  is  also  included  in  Table  1 .  This  loading  mechanism  allowed  the 
spedmen  to  be  tested  in  a  standard  tension  testing  machine  without  the  requirement  of  a 
spedal  loading  fixture  such  as  used  in  Ref.  3.  A  screw-driven  machine  was  used  with  the 
displacement  controlled  at  0.5  mm/min.  Loads  and  displacements  were  digitally  recorded  and 
the  initial  interiaminar  tension  failure  was  generally  recorded  as  a  sudden  decrease  In  load.  To 
ease  the  observation  of  the  location  of  interlaminar  tension  failure  the  sides  of  the  spedmen 
were  painted  white  with  a  water  based  typewriter  correction  fluid. 

Analysis 

To  determine  the  stresses  throughout  the  cross  sedion,  Lekhnitskii  solutions  for  f 
curved  beam  with  an  end  load  and  for  a  curved  beam  with  a  moment  were  used  [6]. 
Corrections  were  made  to  the  Lekhnitskii  solution  to  produce  a  solution  for  a  pure  end  load . 
Corrections  were  also  made  to  the  applied  moment  to  account  for  reductions  in  moment  arm 
length.  An  elementary  beam  theory  solution  and  a  NASTRAN  finite  element  analysis  were 
also  used  to  compare  with  the  Lekhnitskii  solution. 


4 


Lekhnitskil  Solution 

Lekhnitskii  developed  two  sets  of  equations  for  the  stresses  in  a  curved  beam 
segment  with  cylindrical  anisotropy  [6].  A  set  of  equations  was  developed  for  a  cun/ed  beam 
with  an  end  load  and  for  a  curved  beam  with  a  moment  at  each  end.  The  stresses  caused  by 
the  loading  shown  in  Fig.  1  were  calculated  by  superimposing  the  stresses  from  Lekhnitskii's 
two  solutions  (Fig.  3).  However,  the  stress  distribution  applied  to  the  curved  segment  in 
Lekhnitskii's  end  load  solution  produced  an  end  moment  resultant  as  well  as  an  end  force 
resultant.  Consequently,  corrections  were  made  to  the  end  load  solution  by  using  Lekhnitskii's 
moment  solution  to  subtract  out  the  stresses  caused  by  this  additional  moment  (Fig.  3).  Ref.  6 
did  not  indicate  that  an  end  moment  resultant  was  produced  by  the  loading.  Appendix  A 
contains  the  stress  equations  for  the  two  Lekhnitskii  solutions  and  the  modifications  to  the  end 
load  solution. 

Because  Lekhnitskii's  solutions  were  for  a  curved  beam  segment  only,  the  force 
applied  to  the  loading  arm  was  translated  to  the  end  of  the  curved  segment  as  a  moment  and 
a  force  (Fig.  3).  The  stress  equations  are  very  sensitive  to  small  changes  in  the  moment  arm 
length.  Consequently,  corrections  wore  developed  to  obtain  a  more  accurate  moment  arm 
length  [Appendix  B].  Corrections  were  calculated  to  include  the  offset  of  the  loading  pin  from 
the  neutral  axis  (A  +  t!  2,  see  Fig.  1 )  and  the  shortening  of  the  moment  arm  caused  by  the 
rotation  and  displacement  of  the  end  of  the  beam  prior  to  failure.  The  correction  due  to  the 
loading  pin  offset  was  determined  to  be  the  only  significant  effect  relative  to  the  experimental 
error  and  was  the  only  correction  included  in  the  final  moment  calculation.  This  correction 
reduced  the  moment  arm  length  by  up  to  30%  in  some  configurations.  Corrections  due  to  the 
rotation  and  displacement  of  the  end  of  the  beam  were  a  maximum  4.1%.  This  correction  was 
greatest  in  the  thin  specimens. 

The  radial  stress,  including  corrections,  at  any  location  is  given  by  Eq.  1.  The  derivation 
is  contained  in  Appendices  A  and  B. 
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The  analysis  assumed  a  state  of  plane  strain.  However,  the  stresses  calculated  by  assuming 
a  state  of  plane  stress  or  plane  strain  were  almost  identical.  The  variables  in  eq.  (1 )  are 
defined  in  the  appendices.  The  geometry  variables  are  ^so  shown  in  Fig.  1 .  The 
expressions  for  the  tangential  stresses  and  shear  stresses  are  given  in  Appendix  B.  For  the 
loading  arrangement  shown  in  Fig.  1 ,  co  is  equal  to  45°,  and  Eq.  1  reaches  a  maximum  for  any 
value  of  radius,  r,  at  6  equals  45°.  To  determine  the  location,  r,  where  or  is  a  maximum,  Eq.  1 
may  be  differentiated  with  respect  to  rand  equated  to  zero  or  may  be  determined  by 
incrementally  increasing  r.  The  latter  technique  was  employed  in  this  study. 


Elementary  Beam  Theory 

An  elementary  beam  theory  method  to  determine  the  maximum  radial  stress  was 
proposed  in  Ref.  5  as 

L  cos  (co)  -  (a  +  sin(co)j 

2wf(r}ro)i^  (2) 

The  variables  are  defined  the  same  as  those  used  in  Eq  (1 ).  This  expression  accounts  for  the 
maximum  stress  location  being  at  a  position  other  than  the  center  thickness  and  may  be  used 
as  an  approximation  of  the  maximum  radial  stress. 

Finite  Element  Analysis 

An  MSC/NASTRAN  finite  element  analysis  was  conducted  on  the  specimen  to 
determine  the  stress  distribution.  By  using  symmetry,  only  half  of  the  specimen  was  modeled 
(one  loading  arm  and  half  of  the  cxirved  region).  The  force  was  applied  to  the  straight  loading 


arm  so  that  a  comparison  could  be  made  with  the  Lekhnitskii  solution  to  determine  if  the 
translation  of  the  loading  force  to  the  end  cf  the  curved  segment  affected  the  stresses  in  the 
region  of  interest  (6  =  45°).  Four-noded  quadrilateral  elements  were  used  with  24  elements  in 
the  thickness  direction  and  one  element  for  every  degree  in  the  circumferential  direction.  A 
transition  was  made  to  a  coarser  mesh  to  model  the  straight  loading  arm.  For  convenience,  the 
force  was  applied  at  the  neutral  axis  so  a  loading  pin  height  correction  was  not  necessary.  The 
effects  of  geometric  noniinearities  and  large  displacements  were  not  modelled. 

Analytical  Results 

The  radial  and  tangential  stress  distributions,  Cr  and  oe,  in  the  thickness  direction  are 
shown  for  a  3.0  mm  thick  curved  laminate  with  an  inner  radius  of  5.0  mm  at  6  =  45°  and  a 
moment  arm  length,  L,  of  25  mm  in  Fig.  4.  The  shear  stresses  are  zero  at  0  =  45°  and  are  an 
order  of  magnitude  smaller  than  the  radial  stress  for  other  values  of  0.  Hence,  the  shear 
stresses  may  be  considered  negligible.  The  radial  stress  increases  from  zero  at  the  free 
surface  and  reaches  a  maximum  prior  to  the  center  of  the  specimen.  The  tangential  stress  was 
a  maximum  at  the  free  surfaces  and  varied  from  tension  to  compression  as  rwas  increased.  In 
the  region  where  the  radial  stress  is  a  maximum,  the  tangential  stress  is  of  the  same 
magnitude.  Since  the  strength  in  the  circumferential  direction  (fiber  direction)  is  typically  two 
orders  of  magnitude  higher  than  in  the  radial  direction,  a  failure  at  this  location  was  attributed  to 
radial  stress  alone. 

For  comparison,  the  stress  results  from  the  finite  element  analysis  and  from  the 
modified  Lekhnitskii  solution  are  both  shown  in  Fig.  4.  The  maximum  stress  predicted  by 
beam  theory  (Eq.  2)  is  also  shown.  There  were  no  differences  between  the  stress 
predictions  from  the  finite  element  analysis  and  the  Lekhnitskii  solution.  Consequently,  the 
translation  of  forces  in  the  Lekhnitskii  solution  were  assumed  not  to  have  affected  the  stresses 
in  the  region  of  interest  (0  =  45°).  An  analysis  of  a  similar  configuration  indicated  that  the 
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stresses  were  affected  by  the  loading  method  only  within  a  region  that  was  within  1 5°  of  the 
ends  of  the  curved  segment  [3].  The  beam  theory  prediction  was  approximately  8%  lower 
than  the  maximum  stress  predicted  by  the  other  two  methods  and,  hence,  would  give 
conservative  values  of  osc- 

The  distribution  of  radial  stress,  Or,  and  tangential  stress,  oq,  through  the  thickness  for 
beams  of  three  different  thicknesses  is  plotted  at  0  =  45°  in  Fig.  5.  All  the  beams  had  an  inner 
radius  of  5.0  mm  and  a  loading  arm  length  of  25  mm.  The  loads  were  applied  at  the  neutral 
axis  of  the  loading  arm,  and  small  corrections  to  the  moment  arm  were  not  included.  For  a 
given  force,  the  thicker  beams  had  a  lower  value  of  maximum  radial  stress  than  the  thinner 
beams.  Also,  the  location  of  the  maximum  radial  stress  moved  away  from  the  center  and 
closer  to  the  inner  radius  as  thickness  was  increased.  The  maximum  tangential  stress  was 
lower  for  thicker  beams,  and  the  location  where  ae  =  0  did  not  change  significantly  with 
thickness.  Also,  for  the  thicker  48-ply  specimens,  the  curve  is  flatter  near  the  region  of 
maximum  radial  stress  than  for  the  thinner  specimens.  Consequently,  a  larger  percentage  of 
the  cross-sectional  area  is  under  high  stress  for  the  thicker  specimens. 

The  location  of  the  maximum  radial  stress  with  v^ous  thicknesses  and  inner  radii  is 
given  as  a  carpet  plot  in  Fig.  6.  Again,  the  loads  were  considered  to  be  applied  at  the  neutral 
axis  of  the  loading  arm,  and  small  corrections  to  the  moment  arm  were  not  included.  As  the 
inner  radius  increased,  the  maximum  radial  stress  decreased,  and  the  location  of  maximum 
radial  stress  moved  towards  the  center  of  the  thickness.  Additionally,  as  the  laminate  thickness 
increased,  the  maximum  radial  stress  decreased,  and  the  location  moved  away  from  the  center 
of  the  thickn^  towards  the  inner  radius. 

Experimental  Results 

Load-Dispiacement  Curves 

A  plot  of  a  typical  load-displacement  curve  for  each  curve  geometry  is  given  in  Fig.  7. 

In  all  five  of  these  curves,  both  the  loading  arm  and  width  of  the  specimen  are  equal  to 
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25.4  mm.  The  16-ply  specimens  had  an  almost  linear  load-displacement  relationship  to  failure 
with  a  slight  increase  in  stiffness  with  increasing  loads.  This  stiffness  increase  was  probably 
caused  by  a  slight  decrease  in  the  length  of  the  moment  arm  because  of  a  deflection  of  the 
loading  arm.  For  the  24-ply  laminates  with  a  3.2-  or  an  8.5-mm  radius,  the  load-displacement 
curves  were  also  nearly  linear.  For  the  24-ply  specimens  with  the  5.0-mm  radius,  the  stiffness 
was  greatly  reduced  relative  to  the  other  24-piy  specimens  since  the  load  was  applied  using 
hinges  rather  than  the  steel  fixture  as  indicated  in  Table  1 .  These  hinges  were  thin  aluminum 
and  would  bend  on  loading  which  produced  a  nonlinear  load-displacement  curve.  The 
increased  compliance  and  nonlinearity  of  this  configuration  should  not  affect  the  strength 
calculation  since  the  displacements  were  small.  The  48-ply  specimens  were  loaded  using  the 
steel  loading  mechanism  and  were  very  stiff.  Consequently,  the  load-displacement 
relationship  was  nonlinear  because  of  the  small  displacements  relative  to  the  compliance  of  the 
machine  and  loading  fixture.  The  machine  and  fixture  compliance  was  measured  and  found  to 
have  a  similar  nonlinear  behavior  of  the  same  magnitude. 

In  many  of  the  specimens,  subcritical  damage  developed  before  failure  occurred  where 
the  specimen  had  a  significant  loss  of  load  and  stiffness.  The  strength  was  calculated  based 
on  the  load  corresponding  to  this  initial  damage  if  subcritical  damage  developed  prior  to  failure. 
The  development  of  subcritical  damage  caused  small  load  drops  which  resulted  in  steps  in  the 
loading  cunres.  More  of  the  wider  specimens  had  these  steps  than  the  narrower  spedmens. 
An  example  of  a  stepped  loading  curve  is  shown  in  Fig.  7  for  the  24-ply  specimen  with  the 
3.2-mm  radius.  The  cause  for  these  steps  is  discussed  below.  The  1 6-ply  specimens  did 
not  have  any  steps  before  a  sudden  failure  occurred  which  resulted  in  an  approximately  80% 
drop  in  load.  For  the  24-ply  specimens  with  a  radius  of  3.2  mm,  half  of  the  25.4-mm  wide 
specimens  had  steps  while  none  of  the  12.7-mm  wide  specimens  had  steps.  For  the  24-ply 
8.5-mm  radius  specimens,  only  one  of  the  25.4  mm  wide  specimens  had  steps  in  the  loading 
curve,  and  none  of  the  12.7-mm  wide  specimens  had  load  drops.  When  the  24-ply 
specimens  failed,  the  load  dropped  approximately  50%  from  the  maximum.  In  general,  the 
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load  dropped  a  greater  percentage  for  specimens  that  failed  at  higher  loads.  The  48-ply 
specimens  were  unloaded  at  the  first  sign  of  damage  as  indicated  by  a  small  load  drop 
(approximately  70  N)  and  circumferential  crack.  These  damaged  specimens  were  later 
reloaded  until  a  large  drop  in  load  occurred.  Most  of  the  load  curves  had  several  more  steps 
before  a  final  failure  occurred.  The  final  failure  occurred  at  load  that  was  more  than  twice  the 
value  at  which  initial  damage  developed.  At  failure,  the  load  would  typically  drop  by  40%. 

After  failure,  circumferential  cracks  in  the  curved  segment  were  observed  on  the  edges 
of  all  specimens.  These  cracks  extended  slightly  into  the  loading  arm.  Fig.  8  shows  an 
example  of  a  circumferential  crack  for  a  24-ply  specimen  with  an  inner  radius  of  8.5-mm.  The 
16-ply  specimens  tended  to  have  the  most  cracks  with  approximately  seven  cracks  per  side. 
Most  spedmens,  however,  had  one  to  three  cracks  on  each  side.  These  cracks  were  larely  in 
the  same  radial  location  on  opposite  sides  of  the  same  specimen.  The  small  load  steps  that 
occurred  prior  to  failure  resulted  in  a  small  circumferential  crack  on  one  or  both  edges.  For  the 
48-ply  specimens,  the  cracks  formed  by  the  load  steps  always  formed  near  mid  thickness. 

Failure  Strengths 

The  interlaminar  strength,  03c,  was  calculated  for  each  specimen  using  the  modified 
Lekhnitskii  solution  (Eq.  1 )  and  the  load  at  which  initial  damage  was  detected.  The  strengths 
along  with  the  mean  and  coefficient  of  variation  are  given  in  Table  2.  Figure  9  shows  the 
individual  strength  measurement  for  each  of  the  test  spedmens.  The  strengths  were 
expected  to  be  approximately  equal  for  each  data  set.  However,  the  data  sets  generally  fell 
into  a  high  or  a  low  category  of  strength.  The  eight  sets  of  data  that  fell  into  the  low  strength 
category  were  composed  of  all  the  48-ply  specimens  and  the  24-ply  spedmens  with  inner 
radii  of  8.5  mm  and  5.0  mm.  The  high  strength  set  consisted  of  the  16-ply  specimens  and  the 
12.7-mm  wide  24-ply  specimens  with  the  3.2  mm  radius.  The  24-ply  25.4-mm-wide 
spedmens  with  the  3.2  mm  radius  had  the  most  scatter  with  strengths  that  fell  into  both 
categories.  However,  a  significant  amount  of  scatter  occurred  within  each  data  set.  For  each 
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curve  configuration,  the  effect  of  different  width  specimens  was  small.  However,  four  of  the 
five  curve  geometries  showed  a  slight  decrease  in  strength  for  the  wider  specimens.  The 
strength  was  not  significantly  affected  by  the  two  loading  arm  lengths,  L,  used  in  this  study. 

Since  small  sample  sizes  were  used  to  measure  the  interlaminar  strength,  an  estimation 
of  experimental  uncertainty  was  performed  [7].  Many  of  the  variables,  such  as  the  modulus  in 
the  radial  direction,  are  not  known  with  a  high  degree  of  accuracy.  Consequently,  it  was 
desirable  to  know  which  measurements  were  most  crucial  to  the  accuracy  of  the  interlaminar 
strength  prediction.  The  details  of  this  analysis  are  given  in  Appendix  C.  For  the  case  of  an 
inner  radius  of  5.0  mm  with  a  loading  arm  length  of  25.4  mm  and  a  width  of  25.4  mm,  the 
percent  uncertainty  ranged  from  6  6  %  in  the  48-ply  specimens  to  14.0%  in  the  16-ply 
specimens  The  two  variables  that  contributed  the  most  to  the  overall  uncertainty  were  the 
thickness  and  the  loading  arm  length.  The  analysis  was  also  fairly  sensitive  to  the  value  of  the 
inner  radius.  The  overall  uncertainty  was  rather  insensitive  to  material  properties.  A  large 
estimated  uncertainty  of  30%  in  the  radial  modulus, Er.  resulted  in  less  than  a  1 .3%  uncertainty 
in  the  predicted  value  of  radial  stress.  Consequently,  tne  uncertainty  in  the  material  properties 
did  not  affect  the  overall  accuracy  of  the  prediction. 

Fracture  Patterns 

Several  specimens  of  each  configuration  were  cut  perpendicular  to  the  fiber  direction  at 
0  =  45®  using  a  low-speed  diamond  saw.  The  cut  surfaces  were  then  examined  under  a 
microscope  to  observe  the  fracture  pattern  across  the  width  of  the  specimen.  The 
observations  from  the  cut  surfaces  are  listed  below.  In  addition,  the  fiber  volume  fraction  was 
observed  to  vary  through  the  thickness.  This  variation  was  quantified  by  measuring  the  ply 
thickness  variation  across  the  thickness  using  a  microscope.  The  average  ply  thickness  in  a 
small  region  was  determined  by  dividing  a  known  distance  (0.824  mm)  by  the  number  of 
plies  contained  in  that  distance.  The  results  of  these  measurements  are  shown  in  Rg.  10.  In 
general,  the  quality  of  the  laminate  was  highest  in  the  inner  half  of  the  thickness  which 
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corresponded  to  the  region  with  the  smallest  ply  thicknesses.  In  most  laminates,  the  ply 
thicknesses  increased  steadily  from  the  inner  radius  to  the  outer  radius. 

The  16-piy  specimens  with  an  inner  radius  of  5.0  mm  had  the  highest  average  value  of 
interlaminar  strength  (Fig.  9).  These  laminates  also  appeared  to  be  of  the  best  quality.  The 
cross  sections  were  not  observed  to  contain  large  resin  pockets  or  voids  like  many  of  the  other 
configurations.  Also,  the  ply  thickness  did  not  change  significantly  across  the  thickness  of  the 
specimen  (Fig.  10).  The  resin  rich  regions  between  the  plies  were  observed  to  be  very 
small.  The  average  ply  thickness  at  6  =  45°  was  0.123  mm  which  was  the  lowest  of  all  the 
specimen  configurations.  The  fracture  pattern  consisted  of  numerous  cracks  scattered  across 
the  entire  cross  section  (Rg.  1 1 }.  Most  cracks  were  translaminar  and  discontinuous  with  many 
branches  which  gave  the  fracture  pattern  a  "shattered"  appearance.  The  cracking  appeared  to 
be  more  extensive  in  the  higher  strength  specimens  than  in  the  lower  strength  specimens. 
There  was  no  significant  difference  between  the  failure  patterns  in  the  different  width 
specimens. 

The  24-ply  specimens  with  an  inner  radius  of  3.2  mm  also  had  a  high  average  strength 
(Fig.  9).  These  laminates  had  more  than  a  40%  increase  in  ply  thickness  from  the  inner  to  the 
outer  radius.  Consequently,  there  were  large  resin-rich  layers  between  the  plies  toward  the 
outer  radius.  However,  no  large  voids  were  observed.  The  failure  pattern  for  the  25.4-mm 
wide  specimens  consisted  of  several  long  cracks  on  the  outer  half  of  the  thickness  which 
tended  to  be  predominantly  interlaminar.  On  the  inner  half,  the  cracks  were  more  translaminar 
with  many  branches  which  were  similar  to  the  cracks  in  the  16-ply  specimens  but  not  nearly  as 
dense.  No  specific  region  existed  where  the  cracks  were  more  concentrated.  A  lower  and  a 
higher  strength  12.7-mm  wide  specimen  were  examined  and  found  to  contain  fewer  cracks 
than  the  25.4-mm  wide  specimens.  The  higher  strength  specimen  contained  two  interlaminar 
cracks  (interfaces  of  plies  7/8  and  11/12  from  inner  radius)  that  were  connected  by  a 
translaminar  crack  to  create  a  crack  that  spanned  the  entire  width  (Rg.  1 2).  The  lower  strength 
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specimen  contained  a  single  crack  (interface  of  12/13  from  inner  radius)  near  the  edge  which 
was  connected  to  two  cracks  that  ran  two-thirds  of  the  way  across  the  width  of  the  specimen. 

The  set  of  24-ply  specimens  which  had  an  inner  radius  of  5.0  mm  had  an  interlaminar 
strength  that  fell  in  the  lower  group  (Fig.  9).  These  specimens  were  also  used  to  determine 
the  effects  of  a  change  in  moment  arm  length  on  the  interlaminar  strength.  The  ply  thicknesses 
increased  approximately  20%  from  the  inner  radius  to  the  outer  radius  (Fig.  1 0).  On  the  outer 
half  of  the  thickness,  a  large  resin-rich  layer  was  observed  between  every  group  of  four  plies. 
These  large  resin-rich  layers  may  have  developed  during  the  manufacturing  process  where 
the  laminate  was  debulked  every  four  plies.  The  failure  in  ail  of  the  specimens  occurred 
primarily  on  a  single  interface  in  the  resin-rich  layer  between  plies  8  and  9  from  the  outer 
surface  (Fig.  1 3).  The  low  strengths  may  have  resulted  because  of  the  failure  at  this 
manufacturing  defect.  No  difference  in  the  failure  pattern  was  observed  for  the  different 
moment  arm  lengths  or  specimen  widths. 

The  24-ply  specimens  with  an  inner  radius  of  8.5  mm  had  low  strength  values  (Fig.  9). 
These  specimens  also  had  the  highest  average  ply  thickness  (0.1 52  mm)  in  the  curved  region 
of  any  of  the  configurations.  The  laminate  thickness  around  the  radius  increased  by  as  much  as 
25%  going  from  the  edge  to  the  center  of  the  arc  which  resulted  in  a  significantly  higher  resin 
content  than  the  other  24-ply  specimens.  In  addition,  many  small  voids  and  resin  pockets 
were  obsenred  across  the  cut  section.  The  fiber  volume  fraction  also  varied  across  the 
thickness  with  the  thickest  plies  and  largest  resin-rich  regions  in  an  area  two-thirds  of  the  way 
from  the  inner  radius.  Many  interlaminar  and  translaminar  cracks  were  observed  throughout  the 
two  25.4-mm  wide  specimens  that  were  examined.  The  stronger  of  these  two  specimens 
had  more  extensive  cracking,  and  the  largest  cracks  were  located  approximately  45%  of  the 
way  through  the  thickness.  The  weaker  specimen  had  more  voids  and  resin-rich  regions  and 
cfid  not  have  a  predominant  cracking  path  across  the  width.  In  addition  to  these  cracks,  both 
specimens  had  large  translaminar  cracks  that  ran  through  the  resin-rich  region  in  the  outer  half  of 
the  thickness.  The  strongest  of  the  1 2.7-mm  wide  specimens  was  also  sectioned  and 
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examined.  This  specimen  failed  primarily  on  a  single  interface  between  plies  6  and  7  from  the 
inner  radius  and  had  no  cracking  in  the  upper  half  of  the  specimen. 

The  48-ply  spedmens  had  the  lowest  values  of  interlaminar  strength  and  the  least 
strength  variability  (Fig.  9).  The  area  closest  to  the  outer  radius  was  observed  to  contain 
several  large  voids  as  well  as  large  resin  pockets.  The  ply  thickness  increased  1 8%  from  the 
inner  radius  to  the  outer  radius  (Rg  1 0).  The  resin  rich  regions  between  the  plies  were  much 
larger  near  the  outer  radius  than  near  the  inner  radius.  The  individual  plies  could  be  easily 
distinguished  across  the  entire  section  without  polishing  the  surface.  During  loading,  as 
mentioned  above,  all  of  these  specimens  had  many  small  subcritical  failures  which  were 
detected  audibly  and  by  a  small  drop  in  the  load.  Two  specimens  (1 2.7-  and  25.4-mm  wide) 
were  sectioned  after  the  first  sub-critical  failure  was  detected.  In  both  specimens,  three  small 
cracks  were  observed  between  the  ply  interfaces  in  the  mid-section,  and  it  was  these  cracks 
that  caused  the  small  load  drop.  The  photomicrograph  of  the  25.4-mm  wide  section  is  shown 
in  Fig.  1 4a  This  interlaminar  cracking  occurred  primarily  at  the  interface  between  plies  26  and 
27  (from  the  inner  radius)  in  the  12.7-mm  wide  specimen  and  between  plies  28  and  29  in  the 
25.4-mm  specimen.  On  further  loading,  the  cracks  would  accumulate  and  finally  coalesce  as 
shown  by  a  section  of  a  different  specimen  loaded  to  failure  (Fig.  1 4b).  When  the  cracks 
coalesce,  a  more  significant  load  drop  occurred.  Two  specimens  of  each  width  were  sectioned 
following  final  foilure.  The  greatest  concentration  of  cracks  was  located  in  the  inner  half  of  the 
radius.  This  region  contained  numerous  branching  transiaminar  cracks.  The  cracks  in  the  outer 
half  of  the  radius  were  primarily  interlaminar.  Several  cracks  were  also  observed  emanating 
from  large  voids  (Fig.  14b). 


Discussion 

Unlike  the  transverse  tension  specimen  discussed  in  Ref.  8,  the  interlaminar  tension 
specimen  is  not  stressed  uniformly.  The  maximum  radial  stress  occurs  at  a  single  point 
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through  the  thickness  and  around  the  curved  portion.  However,  the  maximum  stress  does 
have  a  constant  value  at  any  point  across  the  width.  Therefore,  any  reduction  of  width  should 
reduce  the  probability  of  large  flaws  in  that  section  and  give  higher  strengths.  In  general,  the 
mean  esc  values  of  the  wider  specimens  were  similar  or  slightly  lower  than  the  narrower 
specimens  indicating  a  small  effect  of  width  on  030.  Also,  the  strength  measurements  were  not 
affected  significantly  by  the  two  loading  arm  lengths,  L,  used  in  this  study. 

Since  the  effect  of  width  and  loading  arm  length  were  considered  insignificant  relative  to 
the  scatter,  a  mean  03c  was  calculated  for  all  specimens  with  the  same  cun/e  geometry  and 
plotted  as  bars  in  Fig.  1 5.  The  two  data  sets  for  the  24-piy  specimens  with  an  inner  radius  of 
3.2  mm  were  not  combined  since  several  of  the  25.4*mm  wide  specimens  had  different  failure 
modes  due  to  the  development  of  subcritical  damage  prior  to  failure.  The  range  of  scatter  is 
also  shown  for  all  the  data  sets.  A  large  range  of  interlaminar  strengths  was  obtained  for  the 
five  curve  geometries.  The  differences  in  the  mean  values  of  interlaminar  strength  are 
indicative  of  laminate  quality  rather  than  geometry  effects.  For  thicker  laminates,  a  larger 
volume  is  under  high  stress,  and  the  probability  of  large  flaws  being  present  is  increased. 
Therefore,  if  ail  the  laminates  were  of  equal  quality  with  an  equal  distribution  of  flaws,  the  thicker 
laminates  should  have  a  lower  value  of  030-  However,  the  effect  of  thickness  cannot  be  readily 
compared  in  the  experimental  results  since  laminate  quality  was  inversely  proportional  to 
thickness.  Consequently,  the  trend  of  decreasing  strengths  with  increasing  thicknesses  must 
be  attributed  to  laminate  quality  which  resulted  in  different  failure  mechanisms.  The  best  quality 
of  laminate  was  produced  for  the  thinnest  specimens  (16-ply),  and  these  had  the  highest 
strengths.  The  48-ply  laminates  were  of  the  poorest  quaiity  and  had  the  lowest  strength.  All 
of  the  48-ply  specimens  had  large  voids  and  resin  pockets  in  the  curved  region  which  were 
visible  to  the  naked  eye.  The  24-ply  specimens  had  a  quality  in  between  that  of  the  1 6-  and 
48-ply  laminates.  The  24-ply  spedmens  with  the  5.0-mm  radius  all  failed  at  the  same  ply 
interfeK»  because  of  a  manufacturing  defect. 
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The  best  indicator  of  strength  for  a  given  laminate  was  the  local  ply  thicknesses.  The 
highest  strength  configurations  had  the  lowest  ply  thicknesses  in  the  inner  half  of  the  thickness 
which  also  corresponds  to  the  region  of  highest  radial  stress.  In  general,  the  quality  of  the 
laminate  decreased  from  the  inner  radius  to  the  outer  radius.  In  Ref  4,  a  similar  result  occurred 
when  significantly  higher  values  of  asc  were  obtained  for  laminates  where  high  compaction 
pressure  was  used  in  manufacturing  curved  spedmens.  The  quality  of  the  laminates  could  not 
be  determined  using  C-scans  since  ultrasonic  equipment  could  not  be  used  in  highly  cun/ed 
regions.  Manufacturing  defects  were  most  likely  the  initiation  sites  for  much  of  the  damage. 

The  strength  appeared  to  be  extremely  sensitive  to  flaws.  The  strength  of  the  lower  strength 
spedmens  was  probably  controlled  by  larger  macroscopic  defects  such  as  voids.  These 
type  of  defects  may  have  been  detected  by  ultrasonic  methods  if  ultrasonic  equipment  could 
be  used  in  highly  curved  regions.  For  specimens  with  higher  strengths,  the  strength  may  have 
been  a  function  of  smaller  microscopic  flaws  such  as  fiber/matrix  debonds. 

Strength  results  from  several  other  studies  were  also  plotted  in  Fig.  15  for  comparison. 
The  mean  value  of  03c  determined  from  flatwise  tension  tests  [2]  fell  between  the  lower  and 
higher  strength  groups.  The  mean  value  of  03c  determined  from  "L"  shaped  curved  laminates 
[3]  also  fell  within  a  similar  range.  However,  the  quality  of  the  laminates  used  in  these  two 
studies  was  not  known  since  no  sectioning  was  done.  Several  values  of  020  [2,3,8]  are  also 
shown  in  Fig.  1 5.  These  values  of  02c  were  higher  than  the  03c  values  measured  in  the  other 
studies  [2,3]  and  in  the  low  strength  group  of  this  study.  However,  these  values  of  G2c  are  still 
significantly  less  than  the  values  of  csc  measured  in  the  higher  strength  group  in  this  study. 

This  result  is  consistent  with  the  results  of  Ref.  4.  The  average  value  of  030  determined  from 
curved  elliptical  specimens  made  of  T300/934  carbon/epoxy  was  almost  twice  the  reported 
value  of  02c  [9]  determined  from  transverse  tensile  tests.  The  curved  elliptical  specimen  was 
manufactured  using  very  high  compaction  pressure  which  produced  a  very  high  quality 
specimen.  T300/934  has  very  similar  properties  to  the  AS4/3501  -6  used  in  the  current 
study. 
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A  Weibull  statistical  analysis  [1 1]  was  conducted  on  a  group  of  higher  strength  data  sets 
and  on  a  group  of  lower  strength  data  sets  [Appendix  D].  The  high  strength  values  of  030 
relative  to  the  reported  values  of  osc  can  be  accounted  for  in  terms  of  volumetric  scaling.  The 
transverse  tensile  tests  involve  a  uniform  stress  over  a  large  volume  of  material  while  the 
interlaminar  tensile  tests  with  cun/ed  beams  have  a  nonuniform  stress  applied  over  a  very 
small  volume.  In  Ref.  8,  Weibull  scaling  laws  were  used  to  scale  data  from  transverse  tensile 
tests  to  specimens  where  only  a  small  volume  of  material  was  under  high  stress.  The  results 
from  Ref.  8  indicate  that  Weibull  scaling  laws  can  accurately  predict  the  difference  in  strength  of 
transverse  tension  coupons  and  the  030  values  of  the  "good"  quality  curved  beams  tested  in 
this  study. 

Strength  data  obtained  from  curved  beam  test  specimens  are  very  relevant  to  the 
design  process.  The  curved  beam  is  a  more  representative  specimen  than  a  flat  laminate  to 
determine  failure  or  delamination  in  a  composite  structure  of  complex  shape.  Specimens  with 
curved  geometries  present  unique  manufacturing  problems  that  do  not  occur  in  fiat  panels. 
Consequently,  failure  modes  and  strengths  obtained  from  curved  beam  specimens  will 
closely  resemble  those  in  the  actual  structure  since  the  test  specimen  will  have  similar 
manufacturing  defects.  Because  of  the  curved  nature  of  the  specimen,  it  is  not  possible  to  use 
conventional  C-scan  methods  to  determine  spedmen  quality.  Hence,  the  only  means  of 
determining  specimen  quality  is  by  testing  and,  as  in  this  work,  obtaining  sets  of  data 
representing  "good"  specimens  and  "bad"  specimens.  Similarly,  in  complex  composite 
structures,  quality  cannot  be  determined  using  conventional  C-scans.  Hence,  all  the  data  from 
the  cun/ed  beam  specimen,  including  the  "bad"  specimen  data,  may  be  useful  in  predicting 
^lure  in  a  structure. 

In  this  study,  interlaminar  tension  strengths  were  highly  dependent  on  the  manufacturing 
quality  of  the  specimens.  This  is  a  severe  problem  since  the  quality  of  highly  cun/ed  parts 
cannot  bo  routinely  determined  using  conventional  C-scans.  Also,  a  high  quality  specimen 
was  difficult  to  manufacture.  However,  severed  alternatives  are  available  to  improve  the 
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through-the-thickness  strength  of  curved  composite  structures.  The  manufacturing  process  can 
be  improved  to  consistently  obtain  a  better  quality  part.  Also,  ultrasonic  methods  can  be 
improved  to  determine  laminate  quality  in  highly  curved  regions.  A  tougher  matrix  may  be 
used  to  reduce  the  flaw  sensitivity  that  a  laminate  with  a  brittle  matrix  has.  A  material  system 
with  through-the-thickness  reinforcement  such  as  textile  preforms  also  may  be  used.  If 
through-the-thickness  reinforcement  is  used,  the  strength  should  be  much  greater  than 
traditional  laminates  and  less  dependent  on  the  manufacturing  process. 

Concluding  Remarks 

A  curved  beam,  loaded  so  that  the  radius  of  curvature  increases,  produces  an 
interlaminar  tension  stress  that  is  zero  at  the  free  surfaces  and  reaches  a  maximum  at  some 
location  through  the  thickness.  The  tangential  stress  is  tensile  at  the  inside  radius  and 
compressive  at  the  outside  radius  and  small  where  the  radial  stress  is  a  maximum.  Therefore, 
if  a  unicfirectional  specimen  is  used  with  the  fiber  direction  in  the  tangential  stress  direction, 
interlaminar  tension  failure  will  oocur.  This  paper  describes  one  configuration  in  which  to  test  the 
curved  beam  as  a  method  to  determine  interlaminar  strength,  030.  The  curved  beam  is  a  more 
representative  specimen  than  a  flat  laminate  to  determine  failure  or  delamination  in  a  composite 
structure  of  complex  shape  since  curved  geometries  present  unique  manufacturing  problems 
that  do  not  occur  in  flat  panels.  The  specimen  contains  a  90*^  bend  with  straight  arms 
extending  from  this  curved  portion.  A  hinged  loading  mechanism  was  attached  to  the  arms 
and  then  pulled  in  tension  until  interlaminar  tension  failure  occurred.  Geometry  effects  of  the 
specimen  were  studied;  namely  the  effects  of  loading  arm  length,  inner  radius,  thickness,  and 
width.  The  specimens  generally  did  not  fail  along  a  single  interface  but  developed  many 
interiaminar  and  transiaminar  cracks.  The  average  strength  from  each  configuration  fell  into  either 
a  high  or  a  low  strength  category,  in  the  lower  strength  category,  failure  may  have  ocojrred 
because  of  large  macroscopic  flaws  such  as  voids  or  resin  pockets  which  were  observed  in 
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failed  specimens.  In  the  higher  strength  category,  the  failure  was  probably  caused  by  smaller 
microscopic  flaws  such  as  fiber/matrix  debonds.  The  030  values  determined  from  the  high 
strength  data  sets  was  significantly  higher  than  the  reported  values  of  020  from  transverse 
tensile  tests.  This  relative  increase  in  strength  can  be  accounted  for  by  using  volumetric  scaling 
laws.  From  the  data  available,  the  loading  arm  length  had  little  effect  on  030.  Width  was  found 
to  only  have  a  slight  effect  on  a^c-  The  wider  specimens  had  slightly  lower  strengths  because 
more  material  is  under  high  stress  and  the  specimen,  therefore,  has  a  larger  probability  of 
containing  a  flaw.  The  effects  of  thickness  and  inner  radius  could  not  be  determined  because  of 
the  large  variation  of  laminate  quality  between  specimens  of  different  curve  geometries.  The 
specimen  geometries  with  the  higher  strengths  tended  to  have  lower  ply  thicknesses  in  the 
inner  half  of  the  thickness  which  corresponds  to  the  region  of  highest  stress. 
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Appendix  A  -  Lekhnitskii  Solution 


Lekhnitskii  [6]  developed  a  solution  for  the  stresses  in  a  curved  beam  with  cylindrical 
anisotropy  for  two  different  loading  cases:  a  curved  beam  under  pure  bending  and  a 
curved  beam  with  an  end  force. 


Case  II  -  Curved  beam  loaded  by  an  end  force  (plane-stress) 

^  - 1  -  pP  ]  sin(0  +  co) 

^  +  ^Hif  ^  *  “I  ■  PP  ] 

''re  pK-^f  - 1  -  pP  ]  COS(e  +  co) 


(A7) 

(A8) 

(A9) 
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where 


-pP)  +  C'  +pP)ln(p) 


(A10) 


(A11) 


The  analysis  was  converted  to  plane  strain  by  replacing  Er,  and  vre  in  Eqs.  A5,  A10, 

^ respectively. 


The  symbols  are  defined  as: 


c  f  Subscript  denoting  modified  stress  equation  for  an  end  force 

Er  Modulus  in  the  radial  direction 

£e  Modulus  in  the  tangential  diiection 
f  Subscript  denoting  stress  induced  by  an  end  force  (Case  II) 

Gre  Shear  modulus 

m  Subscript  denoting  stress  induced  by  a  moment  (Case  I) 

M  Applied  moment 

P  End  force 

r,0  Cylindrical  coordinates  of  any  point  in  the  curved  segment 

/{  Inner  radius  of  curved  segment 

To  Outer  radius  of  curved  segment 

w  Width  of  the  specimen 

Or,  oe,  tre  Stress  components  in  curved  segment 

V  Poisson's  ratio 

0)  Angle  of  the  load  relative  to  the  face  of  the  loaded  end  of  the  curved 

segment 

Dimensions  and  variables  are  illustrated  in  Fig.  1 . 
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Correction  to  Lekhnitskii's  End  Force  Solution 

The  stress  distribution  applied  to  the  end  of  the  curved  segment  results  in  an  end 
force  resultant  and  an  end  moment  resultant.  Consequently,  it  was  necessary  to  derive  a 
correction  to  Lekhnitskii's  end  force  solution  to  subtract  out  the  additional  moment.  This 
correction  is  shown  schematically  in  Fig.  3.  This  correction  is  not  necessary  when  co  equals  to 
zero  as  was  the  case  in  Ref.  3.  This  moment  per  unit  width  (M*)  applied  to  the  end  of  the 
curved  segment  was  calculated  by  integrating,  along  the  face  of  the  segment  end,  the 
product  of  the  tangential  stress  component,  o^,  and  the  distance  from  the  neutral  axis. 


<*  =  [  °  09  (9 . 0)  (/■-  d/-=  sin(co) 

Jn 


(A12) 


The  actual  stresses  in  a  curved  beam  due  to  an  end  force  resultant  only  were  then  calculated 
by  subtracting  the  stresses  caused  by  M*  from  Lekhnitskii's  original  end  force  solution. 
Lekhnitskii's  solution  for  a  beam  under  pure  bending  was  used  to  calculate  the  stresses  due 
to  M*. 


(A13) 

(A14) 

(A15) 

Final  Stress  Equations  for  the  Interlaminar  Tension  Strength  Specimen 

The  stresses  in  the  interlaminar  tension  strength  specimen  were  calculated  by 
superimposing  the  stresses  from  these  two  solutions  (Eqs.  A16-A18). 

Or  =  opi  +  aF«  (A16) 

a0  =  o|‘  +  o|F  (A17) 

t,e  =  T}8  +  x^  (A18) 
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Appendix  B  -  Corrections  to  the  Applied  Moment 


The  stress  analysis  was  very  sensitive  to  the  changes  in  the  length  of  the  moment 
arm.  An  error  in  the  moment  arm  length  translates  directly  into  the  same  percentage  error  in 
the  stress  calculation.  Consequently,  it  was  important  to  calculate  the  exact  length  of  the 
moment  arm  at  failure.  Therefore,  several  corrections  were  considered  for  inclusion  in  the 
basic  moment  calculation.  These  corrections  included  the  offset  of  the  loading  pin  from  the 
neutral  axis  and  the  shortening  of  the  moment  arm  caused  by  the  rotation  and  displacement 
of  the  end  of  the  beam  on  loading. 

The  basic  moment  was  simply  calculated  as  the  force,  P,  multiplied  by  the  length,  L, 
of  the  moment  arm  (B1 ) 

Mbasic  =  p/.sin(45)  (Bl] 


Prior  to  loading,  the  moment  arm  length  was  reduced  by  the  offset  of  the  loading  fixture 
from  the  neutral  axis  of  the  laminate  (Fig.  1 ).  This  correction  was  calculated  as 


M^'*"9»  =  -p(A  +  |)sin(45) 


(B2) 


where  A  is  the  height  of  the  hinge  from  the  surface  of  the  specimen,  and  f  is  the  specimen 
thickness.  This  was  the  most  significant  correction  to  the  basic  moment  calculation.  Omitting 
this  correction  resulted  in  an  error  up  to  30%  in  some  cases.  This  correction  was  more 
significant  in  the  thicker  specimens 

The  next  most  significant  correction  was  in  accounting  for  the  deflection  of  the  end  of 
the  loading  arm  at  failure.  Simple  beam  theory  was  used  to  calculate  the  deflection  caused 
by  the  transverse  component  of  the  load  and  caused  by  the  moment  created  by  the 
loading  point  off^t  from  the  neutral  axis.  A  geometric  nonlinear  interaction  exists  between 
the  tension  load  and  the  transverse  load.  This  interaction  was  accounted  for  using  tabulated 
resK:tion  coefficients  [1 0].  The  tension  load  had  a  negligible  effect  on  the  deflection 
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calcutated  from  the  moment  and  transverse  load  only.  The  deflection  correction  was  the 
most  significant  for  the  16-piy  spedmens  where  an  average  4.1  %  error  would  result  if  the 
deflection  vy^  not  accounted  for.  However,  the  error  was  reduced  for  thicker  specimens 
since  the  moment  of  inertia  increases  proportionally  to  the  cube  of  the  thickness.  The 
average  error  for  the  24-ply  specimens  was  2.5%,  and  the  error  was  less  than  1  %  for  the 
48-ply  specimens.  These  errors  were  considered  insignificant  relative  to  the  experimentai 
error  [Appendix  C]  and  were  not  taken  into  account  in  the  final  stress  equation.  The 
nonlinear  effect  of  shortening  due  to  beam  rotation  was  also  calculated.  This  effect  was 
most  significant  in  the  16-ply  specimens,  but  made  less  than  a  0.2%  in  the  overall  moment 
calculation.  Consequently,  the  effect  of  shortening  was  sdso  neglected. 

The  final  moment  equation  (B3)  that  was  used  was  simply  the  basic  moment  plus 
the  correction  for  the  hinge  offeet . 

/I^total  ^  +  /i^hinge  »  p  j/.  -  A  - 

The  effect  of  end  displacement  and  rotation  only  becomes  significant  for  very  thin 
specimens  with  high  loads  and  does  not  generally  need  to  be  included. 


Final  Stress  Equations  for  the  Interlaminar  Tension  Strength  Specimen 

Eq.  B3  was  substituted  into  Eqs.  A16-A18  to  obtain  the  final  equation  for  the 
stresses  in  the  interlaminar  tension  strength  specimen.  The  final  equations  (B4-B6)  are 
shown  below. 


®r  =  [(-^f  - 1  -  pP  ]  sin(0  +  co) 

p(l  cos  (co)  - 1  (/|  +  To)  sin(co)  -  (a  + sin(co)) 


ro^wg 


1  -  p2K  \  rol  1  .  p2K 


(B4) 

pK+1 
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+  (1  -P)(4fp^-i  -pP]sin(e  +  o)) 

p(t  cos  (m)  -  i-(;i  +  Co)  sin(o))  -  (a  +  |^)sin((0)) 

«* 


(B5) 


1  . 1 


'^re  =  7v^[(“^f  •*■  pK-^f  - 1  -  pP ]  cos(0  +  (O) 


(B6) 


25 


Appendix  C  -  Sensitivity  Study 


An  estimation  of  experimental  uncertainty  was  performed  to  determine  which 
measurements  were  most  crucial  to  the  accuracy  of  the  interlaminar  strength  prediction.  The 
estimation  method  consisted  of  first  estimating  the  uncertainty  of  each  variable  and  then 
applying  a  relationship  that  specifies  how  the  uncertainty  of  individual  variables  propagates 
into  the  overall  uncertainty  in  the  prediction  [7]. 

The  uncertainty  in  each  variable  was  estimated  as  shown  in  Table  C1 . 


TABLE  C1  “Estimation  of  experimental  uncertainties 


Percent  error  for  1 

individual  variables  | 

Hi 

mm 

HIH 

1 

Er 

±  30  % 

0.0 

0.6 

1.3 

2 

Eq 

±  10  % 

0.1 

0.2 

0.5 

3 

Gre 

±  20  % 

0.1 

0.2 

0.5 

4 

Vr0 

±  20  % 

0.0 

0.0 

5 

L 

±  1.00  mm 

3.8 

4.2 

4.7 

6 

P 

±  1  %  (load  cell) 

1.0 

1.0 

1.0 

7 

l\ 

±  0.20  mm 

3.1 

2.8 

2.3 

8 

t 

±  0.20  mm 

13.1 

4.1 

9 

w 

±  0.05  mm 

0.2 

0.2 

maximum  or.  MPa 

87.2 

35.1 

17.1 

overall  uncertainty  (U) 

±12.2 

±3.2 

±1.2 

The  overall  uncertainty,  U,  was  calculated  as 


U= 


(C1) 


where  /  is  the  variable  number,  n  is  the  number  of  variables,  and  5i  is  the  change  in  the 
overall  prediction  due  to  the  uncertainty  in  variable  /.  The  uncertainty  caused  by  each 
variable,  6i,  was  taken  as  the  greater  of 


f(»'i . v\+m . vii)-f(vi . VI . vh)  or  f(vi . vi-uj . Vn)-f(vi . vi . vti) 


where  cr  =  f(vi . vn). 


The  modified  Lekhnitskii  solution  was  used  for  the  experimental  uncertainty  analysis. 
The  uncertainties  in  the  test  results  are  shown  for  three  different  thicknesses,  t,  for  the  case  of 
an  inner  radius,  /),  of  5.0  mm,  a  loading  arm  length,  L,  of  25.4  mm,  and  a  width,  w,  of 
25.4  mm.  The  two  variables  that  contribute  the  most  to  the  overall  uncertainty  are  the 
thickness  and  the  loading  arm  length.  For  a  constant  value  of  the  thickness  measurement 
uncertainty,  the  overall  uncertainty  decreases  with  thickness.  The  overall  uncertainty  caused 
by  the  uncertainty  in  loading  arm  length  decreased  slightly  with  increasing  thickness.  The 
analysis  was  also  fairly  sensitive  to  the  value  of  the  inner  radius,  /].  The  overall  uncertainty 
was  rather  insensitive  to  material  properties.  A  large  estimated  uncertainty  of  30  percent  in 
the  radial  modulus,  E(,  resulted  in  less  than  a  1.3  percent  uncertainty  in  the  predicted  value 
of  radial  stress.  When  terms  are  added  by  the  squares,  only  the  larger  terms  contribute. 
Consequently,  the  uncertainty  in  the  material  properties  did  not  affect  the  overall  accuracy  of 
the  predction. 
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Appendix  D  -  Weibull  Analysis 


A  Weibuti  statistical  analysis  was  conducted  on  a  group  of  higher  strength  data  sets 
and  a  group  of  lower  strength  data  sets.  The  data  from  the  1 6-ply  specimens  (both  widths) 
were  combined  with  that  of  the  24-ply  specimens  with  an  inner  radius  of  3.1 8  mm  and  a 
width  of  12.7  mm  to  form  a  high  strength  group.  A  lower  strength  group  was  formed  by 
combining  the  data  sets  from  the  24-ply  specimens  v\nth  an  inner  radius  of  5.0  mm.  Data 
sets  could  be  combined  since  the  mean  strengths  varied  only  slightly  between  data  sets 
(Table  2),  and  hence  any  volumetric  effects  on  asc  between  the  specimen  configurations 
were  assumed  to  be  small  [8].  Also,  the  specimens  from  the  data  sets  that  were  combined 
were  of  similar  quality  and  had  similar  failure  modes.  The  probability  of  failure  at  a  given 
stress  level,  P(a),  is  given  by 

P(o)  =  1  -  j 

where  m  is  the  Weibull  slope  and  Oc  is  the  characteristic  strength  [1 1].  The  Weibull  slope, 
m,  characterizes  the  scatter  and  may  be  a  constant  for  a  given  material  system.  A  median 
ranking  scheme  was  used  to  determine  the  probability  of  failure  [8].  A  value  of  m  of  7.76 
and  a  value  of  cc  of  84.1  MPa  were  calculated  for  the  high  strength  group.  This  value  of  m 
compares  well  with  an  m  value  of  7.68  obtained  for  1 14  transverse  tension  tests  of 
T300/5208  in  Ref.  1 2.  In  general,  a  large  number  of  tests  are  required  to  determined  an 
accurate  value  of  m  [8].  For  the  low  strength  group,  a  value  of  m  of  4.88  and  a  value  of  Oc 
of  39.1  MPa  were  calculated.  This  smaller  value  of  m  indicates  a  larger  amount  of  scatter. 
This  value  may  be  different  from  the  other  two  values  of  m  since  the  specimens  from  the 
lower  strength  group  all  failed  at  a  macroscopic  manufacturing  defect  rather  than  at  a 
randomly  distributed  flaw  in  the  microstructure  of  the  material. 
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No.  of 
Ries 

Inner 

Radius, 

/l(mm) 

Thickness, 
f  (mm) 

Loading 

Arm 

Length, 

L(mm) 

Rns 

No.  of 
Specimens 

Loading 

Mechanism 

16 

m 

1.96 

25.4 

12.7 

5 

25.4 

6 

fixture 

24 

3.2 

3.35 

25.4 

12.7 

6 

fixture 

25.4 

6 

fixture 

■ 

3.36 

25.4 

12.7 

5 

25.4 

3 

50.8 

12.7 

5 

25.4 

5 

8.5 

3.66 

25.4 

12.7 

6 

fixture 

25.4 

6 

fixture 

48 

B 

6.61 

25.4 

12.7 

6 

fixture 

25.4 

6 

fixture 

Table  1.  Test  Matrix 
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1  No.  of 
1  RIes 

Inner 

Radius, 

/l(mm) 

Loading 

Arm 

Length, 

L  (mm) 

Nominal 

Width, 

>v(mm) 

Interlaminar  Strength,  03c 
(MPa) 

Mean 

CV,  % 

16 

g 

25.4 

12.7 

92.9,  85.5,  66.7, 

73.7.  88.2 

81.4 

13.3 

g 

25.4 

72.6,  77.6.  84.9, 

88.1,  100.4,  62.7 

81.1 

16.2 

3.2 

25.4 

12.7 

89.4,  75.4.  88.1 

65.7,  66.0.  68.8 

75.5 

14.3 

25.4 

18.0.56.2.  76.4. 

39.7.  40.8.  53.4 

47.4 

41.4 

25.4 

12.7 

25.6,  21.6.  32.3.  33.8,  36.6 

30.0 

20.7 

24 

25.4 

36.7,  37.7.  32.8 

35.7 

BEIi 

m 

50.8 

12.7 

37.4.  41.2,34.0.  50.2,  22.3 

37.0 

27.5 

■1 

25.4 

47.6.  34.2,  47.0.  35.7,  39.5 

40.8 

15.3 

8.5 

25.4 

12.7 

48.7,  43.8.  49.9, 

28.8,  44.6,  22.2 

39.7 

28.8 

25.4 

15.9, 14.7,  29.9, 

34.2.  47.8,  32.9 

29.2 

42.6 

48 

g 

25.4 

12.7 

15.8, 16.8, 16.1, 

16.4, 16.7,  18.5 

16.7 

5.5 

■ 

25.4 

17.3, 15.6,  16.8, 

16.9,  18.2,  17.7 

17.1 

5.3 

Table  2.  Interlaminar  Tensile  Strengths 
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Fig.  1 .  >  Schematic  of  interlaminar  tension  strength  specimen 
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Fig.  3  -  Superposition  of  Lekhnitskii  solutions  to  obtain  stresses  for  test  configuration 
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Fig.  4  -  Comparison  of  analytical  data  reduction  techniques 
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-  Effect  of  specimen  thickness  on  stress  cjistribution 


radial  stress  location 
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Maximum  Radial  Stress  Location, 
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-  Typical  load-displaoment  curves 


Fig.  8  -  Circumfeiential  crack  on  the  edge  of  a  24'p!y  specimen  with  an  r/=  8.5  mm 


Individual  Test 


Fig.  9  -  Experimental  interlaminar  tension  strengths 


10  -  Variation  of 
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Fig.  1 1  -  Cross  section  of  interlaminar  tension  failure  of  a  16-ply  specimen 
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Fig.  12  -  Cross  section  of  interlaminar  tension  failure  of  a  24-ply  specimen  with  r  j=  3.2  mm. 
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Fig.  13  -  Cross  section  of  interlaminar  tension  failure  of  a  24-ply  specimen  with  r  j  =  5.0  mm 


initial  damage 
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Fig.  14a  -  Cross  section  of  a  48-ply  specimen  with  initial  damage 
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Fig.  14b  -  Cross  section  of  a  48-ply  specimen  after  final  failure 


Fig.  15  -  Comparison  of  results  with  other  studies 
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